Particle number fractionalization of a one-dimensional atomic Fermi gas with 

synthetic spin-orbit coupling 
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We propose an experimental scheme to simulate the Iractionalization ol particle number by using 
a one-dimensional spin-orbit coupled ultracold lermionic gas. The wanted spin-orbit coupling, a 
kink-like potential, and a conjugation-symmetry-breaking mass term are properly constructed by 
laser-atom interactions, leading to an effective low-energy relativistic Dirac Hamiltonian with a 
topologically nontrivial background field. The designed system supports a localized soliton excitation 
with a fractional particle number that is generally irrational and experimentally tunable, providing a 
direct realization of the celebrated generalized-Su-Schrieffer-Heeger model. In addition, we elaborate 
on how to detect the induced soliton mode with the fractional particle number in the system. 

PACS numbers: 11.27.+d, 67.85.-d, 71.70.Ej 



I. INTRODUCTION 

The idea of fractional particle number (FPN) goes back 
to the Jackiw-Rebbi model [l|, 0] in the relativistic quan- 
tum field theory, where fractionalization of fermion num- 
ber exhibits when a fermionic field is coupled to a topo- 
logically nontrivial background field. The first physical 
demonstration of this remarkable phenomenon was pro- 
posed by Su, Schrieffer, and Heeger (SSH), in which a do- 
main wall in one-dimensional (ID) dimerized polymers, 
such as polyacetylene, induces a zero-energy soliton state 
Q. The particle-hole ambiguity of the zero mode re- 
stricts the fractional fermion number to be only ±i in 
this system [3, Afterwards, achievements have been 
made to generalize it to an irrational fermion number by 
introducing another field to break the conjugation sym- 
metry, such as different on-site energies [6H8|. 

Another famous example of FPN is illustrated in the 
fractional quantum Hall effect regime Q, where the 
Laughin quasiparticles not only have fractional charges 
but also have fractional (anyonic) statistics in two di- 
mensions (2D). Recent search for fractionlization in 2D 
systems has theoretically demonstrated that fraction- 
ally charged excitations may exist in graphenelike [Toj . 
square-lattice [llj and kagome-lattice [12j systems with 
vortex-type order parameters (which describe the mass 
of the analog Dirac fermions in the systems). The newly 
discovered quantum spin Hall insulators were also pro- 
posed for realizing the SSH model based on the proxi mity 
effect, which introduces a magnetic domain- wall [l3l. Il4j. 
Notably, the edge electrons there with the inherent chiral 
symmetry may exhibit a direct signature of FPN [l3| . 

On the other hand, quantum simulation of relativis- 



* Electronic address: 
t Electronic address: 



zwang@hku.hk 



slzhu@scnu.edu.cn 



tic Dirac Hamiltonian by using ultracold atomic gases 
has recently attracted great interest [HI]. For example, 
ultracold fermionic atoms trapped in a honeycomb opti- 
cal lattice (OL) were theoretically proposed to behave as 
massless and massive Dirac fermions 16], and confirmed 
in a recent experiment [T^ |, The atomic gases with the 
synthetic spin-orbit (SO) coupling [l8l - |2~l| through light- 
induced gauge fields (22|, HH were also proposed for in- 
vestigating interesting Dirac dynamics |24l - [29| . These 
cold atom systems provide a highly controllable platform 
for studying a wide range of models in relativistic quan- 
tum mechanics and field theory [15j . Interestingly, Ru- 
ostekoski et al. presented an experimental scheme to re- 
alize [3(| and detect [3l| the fractionalization of particle 
number by using a two-component ultracold Fermi gas 
in a ID optical superlattice. The low-energy effective 
theory for the atoms in the system becomes relativistic 
under certain conditions, and a laser-induced modulation 
of atomic hopping between neighbor lattices with a kink 
profile gives rise to a physical domain-wall, leading to 
soliton modes with FPN. 

Inspired by recent experimental achievements in the 
artificial SO coupling in ultracold bulk bosonic 0, 
and particularly fermionic atoms (20l.[2l|. we here present 
an alternative proposal for realizing the particle num- 
ber fractionalization using a ID atomic Fermi gas with 
the synthetic SO coupling. The required SO interac- 
tions and a kink-like potential are properly constructed 
by dressing atoms with laser beams in the system, such 
that the low-energy fermionic atoms can behave as mass- 
less Dirac fermions coupling to a topologically nontrivial 
background field. As a result, a localized soliton excita- 
tion in the middle of the effective energy gap appears on 
the domain wall, which is a direct quantum simulation 
of the standard SSH model. Another two laser beams 
are used to introduce an effective Zeeman term, which 
shifts the soliton excitation from the zero-energy. For a 
midgap state below the zero energy level, it takes more 
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fractional fermion number from the valance band and less 
from the conduction band, and vice versa for the oppo- 
site case, such that the soliton state exhibits an irrational 
FPN in this general case; moreover, its profile and FPN 
in the system are experimentally tunable. Furthermore, 
we suggest experimentally available methods to detect 
the induced soliton modes with FPN through measuring 
the soliton density distribution and the local density of 
states (LDOS) near the kink. We also discuss the pos- 
sibility of generalizing our proposal to the realization of 
the FPN in higher spatial dimensions. 

The paper is organized as follows. In the next section 
(Sec. II), we propose an experimental scheme to simu- 
late the generalized SSH model with an irrational FPN 
by using a ID SO-coupled atomic Fermi gas. The real- 
ization of a relativistic Dirac Hamiltonian with a kink 
background field is shown, and the induced FPN in the 
system is calculated and explained. In Sec. Ill, we elab- 
orate on how to detect the soliton modes with FPN in 
experiments. Finally, in Sec. IV, we briefly discuss the 
generalization of the system to higher dimensions and 
present conclusions. 




FIG. 1: (Color online) Schematic representation of the laser- 
atom interaction configuration for generating an effective rel- 
ativistic Dirac Hamiltonian. (a) The three ground states |1), 
1 2) and |3) are resonantly coupled to the excite states |e) by 
lasers with Rabi frequencies flj, and two additional lasers Q ly 
and £ls y couples |1) and |3) to |e) with a large detuning Ad. 
The choose for the three ground states is that, |1) and |3) are 
two degenerate Zeeman sublevels which are addressed by laser 
beams with different polarization and |2) is another hypcrfinc 
level with different energy so that it is addressed by a laser 
with a different frequency, (b) The spatial configuration and 
propagating direction of the laser beams. The candidate for 
the fermionic atoms can be 6 Li or 40 K. 



II. SIMULATION OF FPN WITH SO COUPLED 
FERMIONIC ATOMS 

In this section, we show how to simulate the fraction- 
alization of particle number by using atomic Fermi gases 
with the synthetic spin-orbit coupling. Let us start with 
a brief review of the celebrated model describing kink- 
soliton states and arbitrary fractional fermion number in 
the context of relativistic quantum field theory [l|, Q . For 
ID massless Dirac fermions subject to two static bosonic 
fields tpi and tp 2 , the relativistic Dirac Hamiltonian is 
given by 32] 

H D =ca x p x -ip 2 (x)a y + ipi(x)a z , (1) 

where c is the effective speed of light and <J x ,y,z are the 
Pauli matrices. The background field with a kink poten- 
tial is described by @, 0] 

(p 1 (x) = (p1, cp 2 (x -> ±oo) = ±<p%, (2) 

where tp\ and tp\ are constants. The kink ip 2 acts as the 
boundary of two degenerate vacuums Q • The relativistic 
Dirac Hamiltonian with such a topologically nontrivial 
background potential supports an unpaired soliton state, 
which gives rise to fractionalization of particle number 
Q. Moreover, the FPN is generally irrational and takes 
one-half in the standard SSH model with the conjugation- 
symmetry when ip\ is vanishing. 

A. Realizing the relativistic Dirac Hamiltonian in 
cold atom systems 

Now we demonstrate how to realize the wanted Dirac 
Hamiltonian ([1} with a SO-coupled atomic Fermi gas. 



We consider an ensemble of quasi-2D noninteracting 
fermionic atoms with three relevant spin components in 
the ground-state manifold {|1), |2), |3)}, which are reso- 
nantly coupled to a common excited state \e) through the 
standard tripod configuration [33|,[34j as shown in Fig. 1. 
A candidate for the fermionic atoms can be 6 Li or 40 K. 
For 6 Li atoms, the hyperfine levels can be selected as 
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For 40 K atoms, the corresponding hyperfine levels can be 
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The corresponding Rabi frequencies of the three reso- 
nantly coupling laser beams can be parameterized as 

ill = 17 sin a cos8e~ lKX , 

n 2 = ficosae"'^, (5) 

3 = Q, sin a sin 8 e lKX . 

The wave numbers are k and tjk as shown in Fig. 1(a), 
and = \jY^j=i l%| 2 i s tne total Rabi frequency. Here 
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r] = 1 + 5rj with deviation 6r) being for matching the 
resonant-coupling frequency of the second laser beam. 
For the selected atomic hyperfine states in Eqs. (j3|) and 
dlD, 5r] « 5 x 1(T 7 for 6 Li atoms and Srj 3.5 x 1CT 6 
for 40 K atoms. This can be achieved in experiments by 
adjusting the laser frequency. The deviation is negligible 
in our derivations, however, we still use the notation r\ 
in the following discussions for consistency. We further 
adopt uniform plane- wave laser beams that f2, a, and 9 
are all constants, and particularly choose 9 — tt/4. 



respectively determined by 
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FIG. 2: (Color online) Schematic representation of (a) the 
coupling lasers which generate the wanted external potentials 
Vi,2,3 in Eq. <(Sj ; and (b) the spatial configuration of the Rabi 
frequency Ql(x) which forms the needed kink- like potential. 
The frequency difference uii2 = 228 MHz for 6 Li atoms and 
oji2 = 1186 MHz for 40 K atoms for the selected hyperfine 
levels in the text, both of which are much larger than the 
natural linewidth of the excited state |ee) (about 6 MHz). 

The single-particle Hamiltonian for each atom takes 
the form H s = p 2 /2m a + H lnt , where p denotes the mo- 
mentum operator and m a is the atomic mass. The light- 
atom interaction Hamiltonian Hint is given by H m t — 
7i$^j-_i(fij|e)(j| + H.c). Diagonalizing H- mt yields two 
orthogonal dark states 



\Dx) = (e iKX \l) 
\D 2 ) = (c 4M |l) 
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— sina|2). 



(6) 

The dark states are decoupled to the excited state |e) and 
thus are immune to spontaneous emission. They span 
a degenerate subspace, in which the full state of a sin- 
gle atom can be written as \x( r )} = Yj%=i V'»( r )l-Oi( r ))- 
The center-of-mass amplitudes tpi(r) corresponding to 
the spatial wave functions of the two dark states obey 
the Schrodinger equation ihdt$ — H<&, where the two- 
component spinor $(r) = (V>i(r), i/j2(t)) t and the Hamil- 
tonian reads 



and tt„ = <D 3 -(r)|V|A,(r)> with V = £? =1 7,(r)|t)<*| 
[22l l23j | . We consider the adiabatic motion of atoms ini- 
tially prepared in the dark state subspace. It should be 
noted that the two dark states are not the lowest energy 
states in this system so that the adiabatic approximation 
works well only for finite time scales (mainly due to col- 
lisional relaxations), up to several hundred milliseconds 
under realistic conditions [34| . To obtain a lowest-energy 
two-fold (nearly) degenerate subspace, one can adopt the 
optical dressing scheme described in Ref. [13] (where 
atoms with simpler A type configuration) or in Ref. [35| 
(where more atomic internal states and coupling lasers 
are required). 

The Rabi frequencies chosen in Eq.(|5[) can realize the 
first term in Hamiltonian ([T]) with certain potentials; 
however, to implement the required potentials exactly 
in Eq.JT]), more complicated laser configures are needed. 
One possible method to generate the required potentials 
is that we further choose two additional laser beams with 
frequencies v\ and 1/2, as schematically shown in Fig. 
2(a). The first laser beam (denoted by frequency v\) with 
the effective Rabi frequency U L (which takes the real 
Rabi frequency and the detuning into account) is blue 
detuned for atoms in the internal levels |1) and |3), but 
red detuned for atoms in level |2), all of which are far-off- 
resonantly coupled to another excited state |ee). This en- 
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40 K atoms, respectively. In addition, the second laser 
(denoted by frequency 1/2), which is also far-off-resonant, 
is use to create constant energy terms in Vi^ and V2 in 
the following equation ((EJ • The energy difference between 
them [cf. Eq. dHJl] can be realized by detuning the second 
laser from the two-photon resonance with the frequency 
Tin 2 {l — rf cos 2 a)/2m a . Thus the resulting external po- 
tentials are given by 



v 1 = v 3 = m L (x,z) 
v 2 = -nn L (x,z)-^- 



2m„ 



(8) 



2m a 



cos 2 a. 



After introducing all of the laser configurations, we can 
obtain the total resulting potentials for the atoms in the 
laser field as 

A = — fincosa <r x e x + hrjn ( „ ^ ] < - 
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The gauge potential A arises from the position de- 
pendence of the dark states and is given by Aj ; „ = 
i7i(Dj(r)|V|D„(r)) [H, H^. The projecting-induced 
scalar potential and the external potential V are 



V 



Vi 
Vi cos 2 a + V 2 sin 2 a 



with <j) + V = hn^ sin 2 aa z up to an irrelevant constant. 
Note that atoms in such a synthetic non-Abelian gauge 
field behave as electrons with a SO coupling, which can 
be seen from the term p • A in Hamiltonian ([7]). 
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By applying an additional extremely anisotropic trap- 
ping potential to freeze the atomic motions along z axis, 
we arrive at the quasi-lD cases 36]. For ultralow tem- 
perature, the momentum of atoms along x axis p x <C 
%k cos a, such that the p x term in Eq. (J7]) may be safely 
neglected, leading to an effective Dirac Hamiltonian 

H e rj c x a x p x + A{x)a z , (10) 

where c x — fin cos a / m a is the effective speed of light in 
this system and A{x) — TiQl(x) sin 2 a. The Hamiltonian 
(|10p describes a massive Dirac fermion having a position- 
dependent mass A(a;)/c 2 , or in another point of view, a 
massless Dirac fermion coupling to a static background 
field A(x) Q| . If we choose the intensity distribution of 
the laser beam with photon frequency v\ as a kink-type 
function along x axis, then the standard SSH model in 
continuum limit Q is realized in this cold atom system. 
It is interesting to note that the laser-atom interaction 
of A configuration [27j can also be used to realize the 
relativistic Dirac Hamiltonian (|T0|) . and in this case the 
experimental setup can even be simpler. While it is noted 
that such a simplified scheme is unable to be extended to 
realize the generalized SSH model described in Eq. (fTTj) 
below with an irrational FPN. 

To introduce the constant field <p± in Hamiltonian 
([1]), which acts as a mass term and breaks the con- 
jugation symmetry, we can apply two additional laser 
beams to couple the atomic states |1) and |3) to the ex- 
cited state |e) off- resonantly with a large detuning A^ 
as shown in Fig. 1, with the corresponding Rabi fre- 
quencies Qi y = i\n y \e lKX and Cl 3y = \Q y \e~ iKX [Hj]. 
Since | | , | r2 3y | <C Ad, the effective Hamiltonian de- 
scribes a perturbation coupling between states |1) and 
|3), which is given by H p = -ihn p e 2tKX \l)(3\ + H.c. with 
fi p = \rt y \ 2 /Ad 37]. We assume Q p <C SI, so the Hamil- 
tonian H p can not pump the atoms outside of the dark 
state subspace. Mapping H p into the subspace spanned 
by the basis {|-Di), lAa)}, we obtain H p = TiQ p cos aa y . 
Therefore, the total ID effective Hamiltonian for the ul- 
tracold atoms is H\y> — H e + H p . By introducing a uni- 
tary transformation $(x) — > e l ~ (Tx ^(x), we can obtain 
the Dirac Hamiltonian 

Hid = c x a x p x + A(x)a y - Ta z , (11) 

where T = h\£l y \ 2 cos a/Ad. We note that the unitary 
transformation used here is for mathematical convenience 
but involves no manipulation on the system. Compared 
to the original Dirac Hamiltonian (Q}, here the effective 
field — r corresponds to the constant background (p\, and 
the field A(x) should present a kink-like profile, which 
corresponds to — (f2(x). To this end, we can choose the 
spatial profile of the Rabi frequency Q l (x) with the kink 
form as shown in Fig. 2(b), and S1_l(x) = ±S1lo as x — ¥ 
±00. Thus the asymptotic value of A(x) can be denoted 
as A = A(x —> +00) = HQlo sin 2 a. 

Finally in this section, we note that the recent experi- 
ment of generating SO couplings in Fermi gases [20j may 



also be extended to realize the Dirac Hamiltonians (|T0|) 
and ([TT|). In the experiment 20], two spin- 1/2 states are 
chosen as two internal hyperfine states instead of dressed 
states [see Eq. ©] in our scheme, and they are coupled 
by a pair of Raman beams with spatially homogenous 
coupling strength flu. The synthetic SO coupling is just 
one-dimensional with the form p x a z , and there is an ad- 
ditional term related to the Raman co upl i ng SIrct x in the 
single-particle Hamiltonian (see Refs. [18|,[20| for details). 
If one uses Raman beams with spatially imhomogenous 
coupling strength and kink- type profile along x axis [i.e. 
Qr(x) = f2i(ar)], the low-energy effective Hamiltonian 
for the atoms takes the form of the Dirac Hamiltonian 
(ITU1) under a spin rotation. To simulate the generated 
SSH model Hamiltonian (JTTJ) , one needs additional laser 
beams or radio-frequency fields to couple the two spin 
states as a c^-coupling term. 

B. Fractional particle number in this system 

We now turn to calculate the FPN in the proposed sys- 
tem. There is a number of methods for computing the 
FPN of topological solitons Q , including the well-known 
conjugation-symmetry analysis for the zero modes with 
one-half fermion number 0, |38[ . It was first pointed out 
by Goldstone and Wilczek that @, at zero temperature, 
the fractional fermion number of the soliton in this model 
is determined by the kink background field @. The adi- 
abatic condition was imposed there for a valid pertur- 
bation calculation by assuming \dtpi\ <§; m 2 (i = 1,2), 
where m = \/</?i + <~p\ . However, Yamagishi showed that 
the exact result actually does not depend on the adia- 
batic condition [39]. For simplicity, we here employ still 
the adiabatic condition to derive the result in a new but 
simple way. The current of this (1+1)D system is 

f{x) = -{x\tv^G\x+), (12) 

where the Dirac matrices 7 = a z , 7 1 = icr y , 7 s = a x , 
and the Green's function of the relativistic Dirac Hamil- 
tonian ([!]) is given by 



^ l p p - (tpi + i(f2l b ) ' 

Here /i = 0,1 correspond to the time and space compo- 
nents, respectively. 

In the derivation, we have used the adiabatic approxi- 
mation that dipi commutes with p-i^ m -2 and kept the first 
order approximation. Thus the Green's function can be 
written as 

G = [7^ + to - i ^ 2 7 5 )]{G- 1 [7 AI PM + (Vi ~ ^27 5 )]}" 1 

~ [l>*Pv + (<Pl ~ ^27 5 )] (p i} m iyl Ydv{Vl ~ ^27 5 )- 

(14) 
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After a straightforward calculation, given that the chem- 
ical potential (the Fermi level) is zero, we can obtain the 
average current in the background field as Q 



(15) 



where 0(x) = arg(^i + if 2) denotes the angular field of 
the background, and & v is the two- index totally antisym- 
metric tensor. When the chemical potential jl = 0, the 
particle density is given by po(x) — — ^ 9 with re- 
spect to the density of kink-free system. Thus we obtain 
the FPN A/"o = / po(x)dx in this system with <p\ = — V 
and if2 (x) = — A(x) as 



,r 1 f A o\ 1 

Alo = — arctan = arctan 

7r V — r / 7T 



^LoA<j 



sin 2 a 



a, 



■ cos a 



(16) 

It is clear from Eq. (|TB)) that A/"o is generally irra- 
tional and can be an arbitrary fractional number in the 
range (—5, 5) for a finite T. Especially, the conjugation- 
symmetric Jackiw-Rebbi model is obtained in the limit 
r — > 0, i.e., without applying the laser beams VL\ y and 
H,3y. In this case, the soliton is a zero-energy mode with 
one- half fermion number A/o = ±5. It is interesting to 
note that the FPN in this system is widely tunable in 
experiments via laser-atom interactions [cf. Eq. (|16[) ]. 
making it a controllable platform for simulating fraction- 
alization of particle number. 

Fractionalization has been wi dely investigated in rela- 
tivistic quantum field theory @, Ho| and condensed mat- 
ter systems 0, , where it can give rise to interest- 
ing transport phenomena. For example, the existence 
of fractionally charged excitations greatly enhances the 
conductivity in the polymers 0] and may induce quan- 
tized currents in the quantum spin Hall insulators [l3j]. 
The interpretation of FPN in these kink-soliton systems 
is usually in terms of deformation (or polarization) of the 
ground-state vacuum due to the kink which supports a 
single soliton mode [H, |4(| • This fractionalization mecha- 
nism is very different from that in the fractional quantum 
Hall effect regime, where the fractional collective exci- 
tations are described by a strongly correlated Laughlin 
wave function Q. 

To have better understanding of this mechanism, we 
consider another kind of background field with a sim- 
ple but experimentally practical configuration, that is, 
a pair of kink and anti-kink both with a step-function 
profile as shown in Fig. 3(a). Here we assume Ao > 
and a ^> hc x /Ao such that the kink and anti-kink are 
almost decoupled. By solving the energy spectrum of 
Hamiltonian (|11[) at the kink potential (near x = — a) 
with A(x) = AoSgn(a; + a), we find that there is a lo- 
calized midgap eigenstate in the kink at E — T with the 
wavefunction decaying as exp (— Ao|x + a\/hc x ) and the 
energy gap E g = 2y/ Aq + V 2 . It is understood that the 
isolated state picks up a fractional fermion number (i.e., 
FPN) of |7V"o| [see Eq. (JTSJ)] from the effective valance 
band and (1 — |Ao|) from the effective conduction band, 
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t l-\N n 
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FIG. 3: (Color online) Schematic representation of (a) a 
background field with a pair of kink and anti-kink, both 
of which support a localized soliton state with profile ~ 
exp (— Ao|a; ± a\/Thc x ); and (b) the energy spectrum and a 
pair of solitons with energies E = ±F and fractional particle 
numbers ±pVo|- The midgap soliton modes are protected by 
an energy gap E g — + T 2 between the effective con- 

duction band and valence band. At the kink, the soliton state 
picks a FPN of |A/o | from the effective valance band and an- 
other 1 — |7Vo| from the effective conduction band, and vice 
versa for the opposite case at the anti-kink. 



as shown in Fig. 3(b). For an anti-kink potential (near 
x = a) with A(x) = — Aosgn(x — a), the localized soli- 
ton state is obtained at E = —T with the wavefunction 
decaying as exp (— Ao|x — a\/Tic x ) as shown in Fig. 3(a). 
It picks up (1 — \Nq)\ from the valence band and |A/o| 
from the conduction band. For a periodic system, there 
must be pairs of kink and anti-kink. If both states are 
unoccupied, the particle numbers are — at the kink 
and |A/o| — 1 at the anti-kink. When the chemical po- 
tential, i.e., the effective Fermi level in this system, is 
tuned up, the E = — T soliton state is occupied first and 
the particle numbers at kink and anti-kink are =f|A/o|, re- 
spectively. And when both states are occupied, there are 
particles (1 — |Ao|) and |Ao| at the kink and anti-kink. 

From Eq. (fT!)|) . we can see that the FPN Ao depends 
only on the asymptotic value of the kink Ao rather than 
the detailed shape of A(x). In this sense, it is topological 
and is insensitive to local fluctuations of the background 
field. This property enables us to use laser beams of 
different and imperfect intensity distributions compared 
to the ones with the exact kink profile as shown in Fig. 
2(b) and with a wide square-potential profile of a kink- 
anti-kink pair as shown in Fig. 3(a). In experiments, 
the intensity distribution of laser beams can be well de- 
signed and the wanted ones with nearly square-potential 
profiles have been realized [4l|. Although the value of 
Ao is obtained at zero temperature, the corresponding 
FPN for finite temperature T defined as A/r can also be 
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calculated by taking the thermal distribution (i.e., Dirac- 
Fermi distribution) into account |42| . Interestingly, At is 
just depends on the asymptotic value of the background 
field and the temperature T [43 |. At low temperatures, 
i.e., \tpi\P > 1 with /3 = l/k B T (k B is the Boltzmann 
constant), one has At « Ao - e~ lv ° l/3 [42|]. For our pro- 
posed cold atom system with the typical temperature 
T ~ 0.1 /jk and parameter T/h ~ 0.1 MHz, the deviation 



Ao - At 



,-r/3 



-10 



Thus we can conclude that 



the FPN in this system is very robust against the finite- 
temperature modification due to the gap protection. 

Before ending this section, we discuss briefly the mod- 
ifications on the soliton state and its FPN arising from 
the neglected quadratic term in the Dirac Hamiltonian 
(|TU| . For this system of bulk atomic gases, the effec- 
tive cutoff momentum is determined by the Fermi level 
and the typical atomic momentum can be one order less 
than the recoil momentum of photons (for temperature 
T ~ 0.1 jitk and ft cos a ~ 10 7 m _1 ). So we can treat 
it as a perturbation SH — p 2 /2m a . This perturbation 
alters the energy spectrum and also breaks the conjuga- 
tion symmetry, but the Dirac point and the energy gap 
opened by the kink-background remain. From the per- 
turbation calculations, we find that both of the spatial 
profile and the energy of the soliton state are modified. 
For the case of step-function kink potential and when 
r = 0, the spatial wavefunction of the soliton state de- 
cays as [l + (1 + ^4fr^)] ex P (" A o^ + a\/tex), 
which is slightly broader than that in the absence of 5H 
[c.f. Fig. 3(a)]. The corresponding energy is shifted from 
E = to E = ~/^ ) /2m a c 2 . up to the first order perturba- 
tion. As long as this energy shift is very small compared 
to the gap, i.e. Ao/2m a c 2 , <C 1, one can expect that 
the soliton state is robust against the breaking of con- 
jugation symmetry induced by the quadratic term [Toj . 
In this case, the modification of FFN can be estimated 

1 /2m a c 

- - arctan(^ r 

complicated for r ^ 0, but one can still follow a simi- 
lar perturbation argument and obtain the corresponding 
modification of the soliton state at the kink. 



as SAfo 



The situation is more 



III. DETECTION OF SOLITON WITH FPN 

In this section, we propose possible methods for detect- 
ing the fractionalization of particle number in the atomic 
system mainly through the soliton density and the LDOS 
near the kink (or the anti-kink) by using two standard ex- 
perimental detection methods for ultracold atomic gases, 
such as in-situ absorption imaging technique [43[ and spa- 
tially resolved radio-frequency (rf) spectroscopy [iij . 

First, the density distribution of the soliton modes may 
be extracted out from the atomic density measurement 
via optical in-situ absorption imaging [43J. In this con- 
tinuum model, we work in the soliton framework, and the 
physical particle number in the soliton sector is equiva- 
lently defined as being measured relative to the free sector 



without the kink background 0, S3- Thus the density 
distribution of the soliton mode is given by 0, E3 



p (x) = fdE[\y E (x)\ 2 -\ip E (x)f 

= ~f(x)-r(x) 



(17) 



E F =0 



where \& e {^e) and T(x) [T(x)] are the fermion single- 
particle energy eigenstates and atomic density distribu- 
tion in the presence (absence) of kink background f2(x), 
respectively. Note that here we have assumed the ef- 
fective Fermi level at Ep = 0, which can be achieved 
by properly tuning the chemical potential of the atomic 
gas. In this sense, we can measure the spatial density 
distribution of the SO-coupled Fermi gas both with and 
without the kink potential by tuning on and off the laser 
beam f^, which correspond to T(x) and T(x), respec- 
tively. The integration of po(x) gives the value of FPN 
Ao in Eq. (JTHJ) . This detection scheme provides a clear 
physical picture of FPN; however, it is hard to be im- 
plemented in a practical experiment, as there is only one 
atom in the soliton sector (kink and anti-kink) comparing 
to N a — 1 ones in the free sector, where the total num- 
ber of fermions N a is restricted by the chemical potential 
and is usually several orders larger than unit. In addi- 
tion, the number of soliton modes can not be scaled with 
increasing N a . However, the occupation of the soliton 
state affects significantly the atomic density distribution 
near the kink, which may be regarded as a convenient 
feature to identify the existence of solitons. 

An alternative but practical approach to probe the soli- 
ton state is measuring the LDOS p(x, E) near the kink (or 
anti-kink) by using spatially resolved rf spectroscopy (44j , 
which has been proposed to detect other midgap bound 
states in bulk Fer mi g ases 14a , including the zero-energy 
Majorana modes 34. l4(jl The idea is that one first 



uses a probe rf field to induce a single-particle excita- 
tion from the initial state |dj) to an unoccupied fluores- 
cent probe state |ct/), and then imaging the population 
in state \a f) to obtain the spatial information about the 
LDOS [2 lH, H3]. If we assume that the probe field 
is weak and is detuning S r { from the the induced tran- 
sitions, then the population change in state \a,f) can be 
calculated from the linear response theory [45|, l47j 



I(x,E) = — (<zj,(x)a/(x)) cx p ai (x,E- £rf)H(<5rf - E), 

(18) 

where S(-) is a unit step function. For a harmonically 
trapped gas, the chemical potential changes from ft = 
here to jl(x) = ^m a uj 2 x 2 with u> being the trapping fre- 
quency under the local density approximation, and E in 
Eq. (fT5| is thus replaced by E — p(x). Due to the trap- 
ping potential, the energy variation over the length scale 
Iq = y/ E g /m a uJ 2 becomes comparable to the energy gap 
which protects soliton modes. Therefore we could use 
a sufficiently weak trap and experimentally control the 
size of the gap to reach the nonvanishing gap and locally 



7 



homogeneous limit. Thus the previously presented phys- 
ical picture about the soliton mode persists. Since the 
soliton mode in the proposed system has energy T in- 
side the gap and is localized at the kink x = —a, there 
will be a significantly enhanced population transfer with 
frequency 5 r {/h = [V — jl{—a)]/Ti near the kink. The con- 
tribution from the soliton mode would be clearly visible 
and well separated from other quasiparticle contributions 
by the energy gap E g . Thus the soliton density distri- 
bution po(x,T) can be mapped and singled out in this 
way. Compared with the previous detection method, the 
later scheme is insensitive to the fluctuations in the ini- 
tial number of fermions N a since (i) the occupation of 
the soliton mode just depends on the Fermi level (i.e. 
the chemical potential at the kink) and the fluctuations 
of N a will not affect the topology of the Fermi level; (ii) 
the soliton mode is an eigenstate that is robust against 
thermal and local fluctuations in the presence of an en- 
ergy gap. Interestingly, even single atom in |a/) state can 
be detected with the standard quantum jump technique 
if |o/) is selected to be a different hyperfme state [13,13]. 
Furthermore, this rf spectroscopy technique works as an 
analog of the powerful scanning tunneling microscope for 
probing the atomic gases [HI, [45| : another atom will oc- 
cupy the soliton state after the original atom is scat- 
tered out by the probing laser. Therefore, although there 
is only one atom in the kink and anti-kink at a time, 
the population in \a,f) increases with the probing time. 
Therefore, this scheme can be easily implemented in a 
practical experiment. In this method, FPN can be indi- 
rectly deduced from the population of the soliton state 
and then it is actually an indirect method to measure 
FPN. 



IV. DISCUSSION AND CONCLUSION 

Before concluding this paper, we discuss briefly how to 
realize the FPN in 2D and 3D relativistic quantum field 
theories [38| . It has been shown that a 2D Dirac Hamilto- 
nian with a vortex-like spatially inhomogenous mass term 
also supports a zero-energy mode with one-half fermion 
number [hHH H]. The wanted 2D SO coupling act- 
ing as the kinetic term in the Dirac Hamiltonian can be 
generated in the previous laser-atom interaction config- 
uration, such as c x OxVx — c z <j z p z with c 2 = \f)%K sin 2 a. 
Other kinds of SO coupling terms can also be generated 



via appropriate optical dressing [22|. Another crucial 
step is to simulate the position-varying mass term with 
vortex profiles [H| . Specifically in this cold atom system, 
one needs z) with profiles of the vortex type in con- 

trast to the kink type for ID cases shown in Eq. 
Fortunately, the needed laser fields with vortex type de- 
fects can be created by using the so-called optical vortex 
technique [48j , which moreover has been implemented in 
cold atomic gases in experiments [49j . 

FPN can also be present in 3D Dirac systems, where 
the topologically non-trivial background field should be 
replaced by a 3D profile of a magnetic monopole [T], [Hj] . 
For the 3D cases, we need the SO coupling term such 
as cr x Px + CyPy + OzVz , which can be synthesized by the 
atom- light-interaction scheme proposed in Ref. [50j . The 
wanted mass term with monopole profiles may be gener- 
ated by using electromagnetic field superpositions like 
those were used to induce 3D Skyrmions in atomic gases 
& 

In summary, we have proposed an experimental scheme 
to realize the fractionalization of particle number with a 
ID SO-coupled ultracold Fermi gas. A kink-like potential 
and a conjugation-symmetry breaking mass term are con- 
structed by proper laser-atom interactions, leading to an 
effective low-energy relativistic Dirac Hamiltonian with 
a topologically nontrivial background field. As a result, 
a localized soliton mode emerges near the kink, having 
an FPN which is generally irrational and experimentally 
tunable. The proposed cold atomic system provides a di- 
rect quantum simulation of the famous generalized SSH 
model. We have also presented two useful methods to 
detect the induced soliton modes and the FPN in the 
system. In view of the fact that SO-coupled Fermi gases 
were realized in two very recent experiments [20I l2lj , it 
is anticipated that the present proposal will be tested 
experimentally in near future. 
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